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Abstract 

We present one additional study of multi-muon events produced at the Fermilab Tevatron collider 
and recorded by the CDF II detector. We use a data set acquired with a dedicated dimuon trigger 
and corresponding to an integrated luminosity of 3.9 fb"-*^. We investigate the distribution of the 
azimuthal angle between the two trigger muons in events containing at least four additional muon 
candidates to test the compatibility of these events with originating from known QCD processes. 
We find that this distribution is markedly different from what is expected from such QCD processes 
and this observation strongly disfavours the possibility that multi-muon events result from an 
underestimate of the rate of misidentified muons in ordinary QCD events. 

PACS numbers: 13.85.-t, 14.65.Fy, 13.20.Fc 
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This Letter reports on one additional test on the possible origin of multi-muon events 
observed at the Tevatron. These events were identified in a previous study [it of a data 
set acquired with two central (|?7| < 0.7) primary (or trigger) muons, each with transverse 
momentum Pt > 3 GeV/c, and with invariant mass larger than 5 GeV/c^ and smaller than 
80 GeV/c^. That study shows that many long-standing inconsistencies between measured 
and predicted properties of the correlated bb production and semileptonic decay at hadron 
colliders could be explained by the presence of a relevant source of muons which 

appear to be mostly produced beyond the beam pipe of radius 1.5 cm (this contribution is 
whimsically referred to as ghost events because they were unnoticed or ignored by previous 
measurements). Within the large uncertainty of the prediction, mostly based on simulations, 
the observed rate of ghost events is found to be consistent with being produced by muons 
arising from in-fiight-decays of pions and kaons, or punchthrough of hadronic prongs from 
Kg or hyperon decays. However, a search in ghost events for additional muons with > 
2 GeV/c and \ri\ < 1.1 and contained in a cos6 > 0.8 cone around the direction of a 
primary muon selects a small but significant fraction of events with a large content of muon 
candidates that appears difficult to account for in terms of known sources with the present 
understanding of the CDF II detector, trigger, and event reconstruction. 

A more recent study by the CDF collaboration [6] has improved the estimate of the 
contribution of ordinary sources to ghost events. This study addresses in particular the 
contribution from pion and kaon in-fiight-decays. In 1426 pb^^ of data, there are 54437 ± 
14171 ghost events and 12169 ± 1319 ghost events with three or more muons which cannot 
yet be accounted for with ordinary sources. 

In this Letter, we investigate the distribution of the azimuthal angle (50) between the two 
primary muons in events in which both primary muons are accompanied by at least one (or 
two) additional muon candidates in a cos^ > 0.8 cone around their direction, and compare 
it to those for all QCD sources known to produce dimuon events: bb, cc, and T production or 
events in which one trigger muon is due to hadrons misidentified as muons (cosmic rays are 
removed from the data sample and the contribution of secondary interactions in the detector 
volume is negligible jjlj). As discussed in Ref. |l|, known QCD sources produce a handful of 
events with four and none with six muon candidates. However, if the unaccounted multi- 
muon events were generated by a gross underestimate of the number of additional muons 
mimicked by hadrons in ordinary QCD events, the 6(j) distribution of primary muons in 
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multi-muon events would be similar to that of ordinary QCD events in which the large 
contribution of next-to-leading order (NLO) terms due to initial and final state radiation 
results in a broader S(j) distribution than that predicted by the Born (LO) approximation. 
In fact, the 6(f) distribution of pairs of b hadrons or jets is traditionally used to determine the 
relative contribution of NLO to LO terms This type of comparison was also suggested by 
Ref. jsl, in which the excess of multi-muon events is modeled with the decay of two colorless 
particles produced through the exchange of a heavy object. In such a hypothetical case, 
their deviation from the back-to-back configuration in the azimuthal angle {6(f) = it) is only 
caused by initial state radiation of the incoming quarks and is expected to be small. 

The study presented here uses a dimuon data set corresponding to an integrated lumi- 
nosity of 3.9 fb~^ and selected with the same requirements used in Ref. High precision 
charged particle tracking is provided by a large central drift chamber surrounding a trio of 
silicon tracking devices composed of eight layers of silicon microstrip detectors ranging in 
radius from 1.5 to 28 cm in the pseudorapidity region \ri\ < 1 {q]. The tracking detectors 
are inside a 1.4 T solenoid which in turn is surrounded by electromagnetic and hadronic 
calorimeters. Outside the calorimeters, drift chambers in the region |?7| < 1.1 provide muon 
identification. We search events for additional muons using tracks with Pt > 2 GeV/c and 
|?7| < 1.1. The rate of additional muons mimicked by hadronic punchthrough is estimated 
with a probability per track derived by using kaons and pions from D*^ — > tt^D^ with 
— )► K^TT^ decays Il0| . The difference between observed additional muons and 



predicted misidentifications is referred to as real muons. 

The 6(f> distribution for all 3.9 M events is shown in Fig. [H Figure |2] compares to the 
corresponding heavy flavor simulations the 6(f> distribution of trigger muons due to bb and 
cc production. This figure is reproduced from Ref. [lOj] that has measured aij^^i_^^ and 
ac^^^c^^ in a dimuon data set corresponding to a luminosity of 742 pb~^. In the bb case, the 
distribution has an average of 2.5 with a rms deviation of 0.8 rad. The long and important tail 
extending to 50 = is due to NLO terms and the non-perturbative fragmentation function 
of b quarks. In cc events, because of the smaller quark mass, NLO terms are approximately 
a factor of three larger and the fragmentation function is much softer. Accordingly, the 6(f) 
distribution has a smaller average (2.4 rad) and a larger rms deviation (0.9 rad). 

The azimuthal-angle distribution for primary muons produced by T(IS') decays is ex- 
pected to be similar to those for heavy flavors because the final state contains a bleaching 
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FIG. 1: Distribution of the azimuthal angle S(j) between the two trigger muons for all events. 



gluon recoiling against the T meson. This distribution, shown in Fig. |3l is constructed using 
muon pairs with invariant mass in the range 9.28 — 9.6 GeV/c^. As in Ref. 10|], the combina- 
torial background under the T(IS') signal is removed with a sideband subtraction technique. 
A similar 6(f) distribution is also expected for those cases in which one muon is mimicked by 
a track in the jet recoiling against a muon due to a heavy-quark semileptonic decay. Figure [3] 
shows the 6(j) distribution of primary muons when one of them is mimicked by pions pro- 
duced by Kg decays. As in Ref. {g], we select Kg — )■ tt+tt^ with a. -k ^ jj, misidentification 
by combining primary muons with tracks of opposite charge and px > 0.5 GeV/c. We select 
pairs consistent to those arising from a common three-dimensional vertex. We also take 
advantage of the K^ long lifetime to suppress the combinatorial background. We further 
require that the distance between the K^g vertex and the event primary vertex, corrected by 
the K'g Lorentz boost, corresponds to ct > 0.1 cm. We select K^ candidates with invariant 
mass in the range 0.47 — 0.52 GeV/c^ (see Fig. 3 of Ref. j^), and remove the combinatorial 
background with a sideband subtraction technique. 

In summary, the 5(j) distributions of primary muons produced by known QCD processes 
peak at 5(f) ^ vr, and exhibit a significant tail extending to 5(j) = 0. Depending on the 
production mechanism, the mean and rms deviation of these distributions are in the range 
of 2.4 — 2.5 rad and 0.7 — 0.9 rad, respectively. 
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FIG. 2: The distributions (•) of the azimuthal angle 6(p between trigger muons due to (left) bb and 
(right) cc production are compared to the corresponding heavy flavor simulations (o). Distributions 
are normalized to unit area. 



The 6(f) distributions in the subset of events in which each trigger muon is accompanied by 
at least one or at least two additional real muons are shown in Fig. HI These 6(f) distributions, 
with mean of 2.9 rad and rms deviation of 0.2 rad and without any tail below 6(f> = 2.5 rad, 
are different from those of primary muons due to all known QCD sources. 

In conclusion, as mentioned earlier, within our present understanding of the CDF-detector 
response no known sources produce events in which each cos^ > 0.8 angular cone around a 
primary muon contain at least two additional real muons. Had the additional muons been 
produced by a subtle failure of our method to evaluate the fake-muon contribution, the 
resulting 6(f> distribution of primary muons would have been found consistent with those 
typical of ordinary QCD processes. 
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removed with a sideband subtraction method. 
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